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ABSTRACT 



/Kms. We show how the build-up of magnetic gradients in the Sun's corona may be inferred directly from photospheric velocity 
data. This enables computation of magnetic connectivity measures such as the squashing factor without recourse to magnetic field 
extrapolation. 

Methods. Assuming an ideal evolution in the corona, and an initially uniform magnetic field, the subsequent field line mapping is 
computed by integrating trajectories of the (time-dependent) horizontal photospheric velocity field. The method is applied to a 12 
hour high-resolution sequence of photospheric flows derived from Hinode/SOT magnetograms. 

Results. We find the generation of a network of quasi- separatrix layers in the magnetic field, which correspond to Lagrangian co- 
herent structures in the photospheric velocity. The visual pattern of these structures arises primarily from the diverging part of the 
photospheric flow, hiding the eff'ect of the rotational flow component: this is demonstrated by a simple analytical model of photo- 
spheric convection. We separate the diverging and rotational components from the observed flow and show qualitative agreement with 
purely diverging and rotational models respectively. Increasing the flow speeds in the model suggests that our observational results 
are likely to give a lower bound for the rate at which magnetic gradients are built up by real photospheric flows. Finally, we construct 
a hypothetical magnetic field with the inferred topology, that can be used for future investigations of reconnection and energy release. 

Key words. Magnetic fields - Sun: photosphere - Sun: corona - Sun: magnetic topology 



1. Introduction 

This paper proposes a straightforward method to study the time- 
dependent build up of structure in the Sun's coronal magnetic 
field, based on observations of horizontal velocity fields in 
the solar photosphere. The ultimate objective is to determine 
w^hether the energy built up and released in the coronal magnetic 
field as a result of photospheric convec tion i s suflftcient to heat 
the corona via the Parker mechanism (Parker 1972). In this the- 
ory, footpoint braiding generates localised magnetic gradients 
and thin current sheets in the corona, leading to ubiquitous re- 
connection. There are two requirements that must be satisfied 
by the observed motions if this is to work: (1) they must have a 
tendency to create the required magnetic gradients, and (2) they 
must do so quickly enough to generate a sufliicient overall recon- 
nection rate. The method proposed here aims to investigate these 
questions in a practical way, given currently available observa- 
tions. 

Numerous studies have found that the topology and connec- 
tivity of the coronal magnetic field play a primary role in deter- 
mining when and w here magnetic reconnection will take place 
(iBirn & Priestll2QQ7h . An increasingly popular means to charac- 
terise the connectivity of 3D magnetic fields has been through 
the s o-called squa shing factor Q of the magnetic field line map- 
ping (^Titov et an20Q 2). Quasi- separatrix layers (QSLs),where Q 
is high, represent locations in the magnetic field with large gradi- 
ents in field line connectivity. They have been identified with the 
locations of flares and X-ray bright points ( Demoulin et al. 1996, 



11997b iMandrini et al.|[T996t iGaizauskas et al.lll998l: fWang et al.l 
'2000^. Magnetic reconnection is suggested to occur preferen- 
tially at QSLs (Demoulin 2006; Santos et al. 2008). 

Although Q is defined solely by the connectivity of field lines 
(the mapping between photospheric footpoints), existing studies 
have all calculated Q by first constructing a 3D magnetic field, 
then tracing field lines to determine the mapping. Our method 
determines the field line mapping, and hence Q, in a fundamen- 
tally diflTerent way. Given some initial field line mapping at time 
to, and assuming an ideal evolution of the coronal magnetic field, 
the field line mapping at a later time depends only on the se- 
quence of photospheric footpoint motions. We propose to take 
advantage of this fact to compute the field line mapping directly 
from photospheric velocity data. 

The main advantage of our proposed method is that it avoids 
the need to extrapolate a 3D magnetic field from photospheric 
magnetograms. This is problematic for studies of magnetic 
topology because standard techniques (such as potential fields) 
applied to a sequence of photospheric magnetograms will not 
give the correct field line topology commensurate with an ideal 
evolution from one time to the next. To avoid this problem, a 
number of studies have used time-dependent 3D simulations to 
model how topologi cal structure develops in response to simple 
boundary motions (iMilano et al.1 [1999b iGalsgaard et al.l [20031: 
lAulanier et al.l2005l:lMasson et al.l2009l) or to bou ndary motions 
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derived directly from observed magnetograms (Mackay et al.l 
2011). But such time-dependent simulations inevitably suflfer 
from numerical dissipation, leading to inaccuracies in magnetic 
topology. Our method avoids this. 



1 



A. R. Yeates et aL: Photospheric flows and coronal magnetic structure 



The main limitation of our proposed method is that it can- 
not determine the initial field line mapping at time to, only 
that resulting from subsequent footpoint motions. The initial 
mapping will likely contain pre-existing magnetic structure and 
QSLs. Unfortunately, existing extrapolation techniques cannot 
uniquely determine this initial magnetic topology from available 
magnetogram observations. Note that this limitation is shared by 
all of the studies cited in the previous paragraph, which must also 
assume some initial magnetic field; typically a potential field 
extrapolation is used. A similar limitation applies to studies of 
magnetic helicity using p hotospheric flows and magnetograms 
("Pemo ulin & Pariaj 12009). In this paper, we simply assume a 
uniform initial field, whose field line mapping is the identity. 
This is the most conservative choice: Q values derived from the 
subsequent mappings will almost certainly be a lower bound for 
the real coronal magnetic field. 

An advantage of taking the initial field to be uniform is that 
one may construct a hypothetical magnetic field with the cor- 
rect topology simply by taking the magnetic field lines to be the 
trajectories of the velocity field. Examples of such trajectories 
are shown in Fig. [H where time increases vertically. The cross- 
sections show the magnitude of for this hypothetical magnetic 
field; these may be thought of as the photospheric magnetograms 
at subsequent times. (See Appendix lAl for details of the calcula- 
tion.) The concentration of magnetic flux in a "network" of con- 
vective cell boundaries is evident. Of course, this hypothetical 
field does not represent a realistic extrapolation of the magnetic 
field that would be seen on the Sun, since it starts from a uniform 
field at to, and only takes into account photospheric motions at 
one end of the field lines. Rather, we envisage using it as the 
starting point for 3D MHD simulations investigating energy re- 
lease, with the advantage of having determined accurately the 
change in field line mapping. 

To demonstrate the proposed method, we apply it to a 
12 hour sequence of photospheric velocities derived by local 
correlation tracking in Hinode/SOT (Solar Optical Telescope, 
iTsuneta et al.ll2008b magnetograms. The observed velocities are 
described in Sect. [2l while the inferred magnetic field line map- 
ping is presented in Sect. O In addition to Q, we compute the 
finite-time Lyapunov exponent (FTLE) field cr. This measure is 
a popular method in fluid mechanics for identifying so-called 
Lagrangian Coherent Structures (LCSs) in velocity fields. Like 
2, cr measures the maximum separation rate of initially nearby 
trajectories, and we illustrate how QSLs in the field line mapping 
correspond to LCSs in the photospheric velocity field. In Sect.|4] 
we explain the pattern observed in the 2 or cr fields using a sim- 
ple analytical model of photospheric convection. By varying the 
model parameters, we predict how the field line mapping would 
be expected to change given observations at higher resolution of 
faster flows. Conclusions are given in Sect. [3 

2, Photospheric Velocity Data 

Our velocity data have been derived by local correlation track- 
ing in magnetograms, although the method could be applied to 
velocity fields from any sourc e. Detailed analysis o f the data re- 
duction procedure is given by I Welsch et al.l (1201 ll) . Briefly, we 
use Stokes V// from Hinode/NFI (Narrowband Filter Imager) 
observations in Fe I 6302A of active region 10930. These wer e 
calibrated to gauss following Equation (1) of llsobe et al.l (l2007l) . 
and the noise level estimat ed at ~ 17 G b y fitting the core of his- 
togrammed field strengths (iHagenaar et al. 1999). In view of the 
subsequent reduction of noise by averaging in the tracking pro- 
cedure, a tracking threshold of 15 G was chosen, with no veloci- 




Fig. 1. A hypothetical magnetic field (Appendix lAb whose field 
line mapping matches that inferred from the observed velocity 
field. The magnetic field lines in this particular field are simply 
trajectories of tracer particles in the 2D, time-dependent velocity 
field (with time increasing vertically). Contour slices show at 
diflferent z with the colour scale in gauss: the initial distribution 
(on the lower boundary) is a uniform field 5^(x,j, 0) = 88 G, 
which was the average field strength in this region of the original 
SOT magnetogram. 

ties assigned to pixels below this threshold. The magnetogram 
pixels are binned (2x2) from 0.16'' to 0.32'', consistent with 
SOT's 0.3" diflfraction limit at this wavelength. The cadence of 
the images is ~ 121s, and the sequence runs from 14:00UT on 
12 December to 02:58UT on 13 December 2006. 

The velocity field is extracted from the magnetograms 
u sing the Fourier l ocal correlation trac king (FLCT) method 
( Welsch et al."'2004; Fisher & Welsch 2008). The method has a 
number of parameters: optimum values have been determined by 
an autocorrelation analysis, aiming to maximise frame-to-frame 
correlations a nd ensure robustness in the velocity estimate (see 
I Welsch et al.i r2011). Here, the windowing/apodization parame- 
ter is set to 4 to avoid too much spatial averaging of small-scale 
flows. The sampling time between subsequent frames is cho- 
sen as = 8 mins. This is small enough to avoid significant 
decorrelation, but large enough to allow for boxcar averaging of 
5 magnetograms to produce each frame, which greatly reduces 
noise. We have repeated the calculations with A^ = 4 mins with 
qualitatively similar results. 

For the analysis in this paper, we select a unipolar plage re- 
gion of size 12.4 Mm x 12.4 Mm (approximately 17" x 17"), 
away from the main sunspots, as shown in Fig. [2] (left). This is 
to avoid the large-scale flow associated with emerging flux and 
rotation of the sunspots. Since the magnetic flux in our region 
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is concentrated in the supergranular lanes, there are inevitably 
areas where the line-of- sight magnetic field is too weak for re- 
liable estimation of the velocity. This particular region has been 
chosen to minimise this problem over the length of the time se- 
quence, although there are several regions where the velocity 
suff'ers locally from high-frequency noise. We have removed this 
noise with minimal disturbance to the well-resolved regions by 
applying a low-pass (Butterworth) filter to the velocity fields in 
Fourier space. Histograms of the velocities both with and with- 
out filtering are shown in Fig. [2] (right). The mean flow speed 
is of the order 0.1 kms"\ which is rather lower than reporte d 
speeds for granular flows (~ 1 km s~^ Rieutord & Rinconl2 010V 
There are a number of possible reasons for this. Firstly, there 
is a likely averaging eflTect due to the convective cells being 
close to our spatial resolution of 0.3''. In addition, comparative 
tests show that FLC T has a bias toward underestimating speeds 
(IWelsch et al.ll2007l) . However, it should be noted that the FLCT 
method tracks coherent magnetic features, which are expected to 
move more sl owly than surrounding plasma due t o suppression 
of convection ("Title e t al.lll992t 'Berger et al.lll998h . The possible 
eflTect of faster flows is explored in Sect. IH 

Figure [3] shows snapshots of the horizontal divergence of the 
(filtered) velocity. 



12-Dec-2006 14:04 



A = V • V = dvjdx + dvyldy, 
and vertical component of the vorticity, 

cl) = • (V X v) = dvyldx - dvxidy. 



(1) 



(2) 



The root-mean- square (RMS) divergence over the whole 12 hour 
dataset is 2.06 x lO'^s'^ (or 3.49 x lO'^s'^ before filtering), 
while the RMS a; is 1.97 x 10"^ s'^ (3.16 x 10"^ s'^ before filter- 
ing). Thus the divergence and curl of the velocity field are com- 
parable in magnitude. This will be important for understanding 
the inferred magnetic field line mapping. 



3. Inferred Magnetic Field Line Mapping 

The field line mapping is simply given by following trajecto- 
ries (particle paths) in the observed time-dependent 2D velocity 
field. So a particle starting from (xq, jo) at time is mapped to 
f(xo,jo,^i) = (/x(-^o,Jo,^i),/);(-^o,Jo,^i)) at time where f is 
found by integrating 



dt 



= v(f(xo,jo,0) 



(3) 



from t = to to t = ti, for which we use a second-order method. 
The observe d velocity fields are interp olated using a local tricu- 
bic method ('Lekien & Marsdenll2005l) that also gives continu- 
ous first derivatives of v. Linear interpolation does not produce 
smooth enough trajectories for computing the Lagrangian struc- 
tures we are interested in. Equation ([3]) is integrated on a grid 
of 864 X 864 starting points to give the 2D field line mapping at 
a given end-time ti . The high resolution is needed to accurately 
determine the Q and cr fields, which typically vary on a smaller 
scale than the velocity field itself (see lShaddenll201 ID . 

To analyse the resulting time sequence of mappings for dif- 
ferent ti, we compute two measures of the mapping gradient: 
the squashing factor Q and the FTLE field cr. The former is fre- 
quently used to characterise 3D magnetic field structure, while 
the latter is used to characterise particle paths in time-dependent 
2D velocity fields. In fact, both are rather similar measures of 
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Fig. 2. Left: Location of the 12.4 Mm x 12.4 Mm analysis region 
(black box) in the full SOT line-of-sight magnetogram (V = 
725 km). Right: Histogram of measured |v| in the analysis region, 
over the entire 12 hour dataset. The solid line shows the filtered 
data, while the dashed line shows the original unfiltered data. 
Each curve is normalised by its own maximum. 



the local rate of stretching at a given point, and both are defined 
in terms of the Jacobian matrix 



J(xo,yo,ti) 



'dfx/dxo dfx/dyo\ ^((^ b 
dfy/dxo dfy/dyoj \c d 



(4) 



3.1. Quasi-separatrix Layers 

The squashing factor at a point (xq, yo) of the mapping from time 
to to ti is 

T 

Q(xo,yoJi)= (5) 

where T = Tr(/^/) = a^+b^+c^+d^ andZ) = det(/) = ad-bc. It 
is a dimensionless measure of the local degree of stretching and 
squashing of a magnetic flux tube under the field line mapping 
(iTitov et al .1120021) . Large values of Q signify the locations of 
strong gradients in the mapping. We find computation of Q to be 
more robust using the equivalent form 



^ R 



(6) 
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Divergence 



Curl 




Fig. 3. Snapshots of A (left) and oj (right) in the analysis region, 
at 14:04UT. The grey scale extends to the maximum absolute 
value of each quantity. 



where R = ^JA+/A- is the ratio of singular values of / (i.e., 
A+ are the eigenvalues of the positive-d efinite matrix J^J; 
iTitov e t al. 2002; Richards on & Finnll201lh . 

We have computed 2 as a function of ti using the inferred 
field line mapping, and the resulting distribution of log^g Q is 
shown in Fig. |4]at two diff'erent times. The logarithmic scaling 
is introduced because certain trajectories typically become ex- 
ponentially separated in time: these are precisely the locations 
of strong gradients in the resulting field line mapping. These 
thin "ridges" of high Q, usually known as quasi- separatrix lay- 
ers (QSLs), are visible at a number of locations in the 12-hour 
snapshot. They are interspersed with more diff'use regions of Q. 
As a simple measure of the overall structure, we show the time 
variation of J log^g Qdxdy as the thick, solid curve in Fig. (6] 
This shows an approximately linear increase through the twelve 
hours of observations. 

All of our figures are shown in the "initia l" frame (xp, yp), so 
that Q SLs correspond to unstable manifolds (iRichardson & Find 
I2OIII) . By definition, Q is independent of the direction of map- 
ping along a particular field line, but when calculated for the 
inverse mapping and plotted at (fx(xo,yo,ti), fy(xo,yo,ti)), the 
pattern would diff'er and the apparent QSLs would then corre- 
spond to stable manifolds. These stable manifolds tend to be 
visible observationally because they form the "network" along 
which magnetic flux concentrates over time. This eflTect is visi- 
ble in Figure [U Conversely, we would expect an anti-correlation 
between the pattern of photospheric magnetic flux at time ti and 
the locations of QSLs in Figure IH which correspond to unstable 
manifolds. 



3.2. Lagrangian Coherent Structures 

An alternative measure of the mapping gradie nt is given by 
the finite-time Lyapunov exponen t or FTLE (iHallei] 1200 U 
IShadden et al.ll2005l: IShaddenll201 Ih . defined as 



cr(xo,yoJi) = 



In 

ki -^oi' 



(7) 



where as before A+ is the largest eigenvalue of J^J. The quantity 
cr gives the separation rate between two initially nearby trajec- 
tories: if two points are initially separated by a small distance 
1^0 1 at time to, then their separation at a later time ti will be 
^ Qxp[cr(ti - ^p)]|fp|. We shall use only the forward time 
FTLE, although one may define also a backward FTLE using 
the inverse mapping. 

In fact, comparing the definition of cr to that of Q in Equation 
([5]) shows that both quantities measure essentially the norm of 




4 6 8 
Time (hrs) 

Fig. 6. Integrated squashing factor, J log^Q Qdxdy as a function 
of time for the observed velocity field (thick lines) and analytical 
models of Sect. |4l (thin lines). The solid lines correspond to the 
observed v and combined model 3. Dotted and dashed lines cor- 
respond respectively to Vdiv (or model 1) and Vrot (or model 2). 
The given slopes are for a least-squares linear fit to the observed 
data. 

the matrix J (which is the Cauchy-Green tensor of the de- 
formation imposed by the field-line mapping). The quantity Q 
uses the Frobenius norm of the matrix, whereas cr uses the spec- 
tral norm. They differ further in that cr includes the logarithm 
in its definition, and has units of inverse time, whereas Q is di- 
mensionless. In the particular case that D = 1, it follows, in the 
high- 2 layers, that Q ^ \A+\, so that we may make the direct 
correspondence InQ ^ 2\ti - to\cr. 

The distribution of cr is shown in Fig. [S] for the same two 
times ti as Fig.lH We see the development of thin ridges of large 
(positive) cr, which represent the locally strongest repelling ma- 
terial surfaces in the flow. These are interspersed by regions of 
negative cr, indicating converging trajectories. The ridges in Fig. 
[5] coincide precisely with the quasi- separatrix layers in Fig. lH 
in the fluid dynamical context they are known as Lag rangian 
Coherent Structures (Haller & Yuanll2000l: IShadden et al.ll2005D . 
They are termed "lagrangian" because they are defined by the 
fluid motion rather than an instantaneous snapshot, and "coher- 
ent" because they have distinguished stability compared to other 
nearby material surfaces. Analysis of such structures has become 
an important tool in fluid dynamics because they often reveal 
the mechanisms underlying transport in complex fluid flows: 
for example, the patterns traced out by visual markers such as 
dye. With our assumption of ideal MHD, the coronal magnetic 
field inherits these topological structures directly from the pho- 
tospheric flow. 



4. Interpretation 

The spatial patterns of Q and cr in Figs. |4] and [5] reveal the 
topological structure of the coronal magnetic field generated 
by the observed photospheric velocities. We can understand 
the origin of this structure using a simple ana lytical model of 
two-dimensional convection, similar to that of ISimon & WeissI 
(1989). We will demonstrate the eff'ects of both diverging and ro- 
tational flows, before using the model to predict how the inferred 
field line mapping might diff'er using higher-resolution velocity 
observations. 



4.1. Analytical Model 

We present three models, each computed on a grid of size 
12Mmx 12 Mm. 
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Fig. 4. Contours of log^o Q ^^^^^ twelve 
hours (at 20:24UT and 02:24UT respectively), 
with initial time to at 14:04UT (available as a 
movie). 



> 



12-Dec-2006 20:24 



13-Dec-2006 02:24 



6 8 
x(Mm) 



E 

4^ 6 



■ ^ 1 

\ 




■ 

yf 


■ 

- 



I 



0.10 
0.07 

0.03^ 
to 

b 

-0.00::^ 

LU 

_l 

-0.03 1 

-0.07 
-0.10 



X (Mm) 



Fig. 5. Finite-time Lypanunov exponent 
(FTLE) field cr at times 20:24UT and 02:24UT 
(available as a movie). 



1 . Diverging flow. Here the velocity field v is a superposition 
of convective "plumes" 



(x-Xt) ^ (y-Yt) 



Vri, (8) 



where (X/, 7/) is the centre of the /* plume, = 
and the plumes have a purely radial 



^(x-Xt)^^(y-Yt)\ 
velocity with respect to their centre of 



(9) 



Thus a; = and A 9^ 0. The centres (X,, 7,) of the = 81 
plumes are selected at random from a uniform distribution, 
over a larger domain -3 < Xi < 15,-3 < Yi < 15 to 
avoid boundary eflTects. To approximate the observations, 
the peak velocities Vi are selected from a normal distribu- 
tion with mean (V/) = 0.1 km ^"^ and standard deviation 
sd(Vi) = 0.03 km ^"^ while the widths Ri are selected from 
a normal distribution (Ri) = 1 Mm and sd(Ri) = 03 Mm. 
For computational convenience, the same pattern of plumes 
is used for a "coherence time" of 15 minutes, before a new 
pattern is chosen. (Qualitatively similar results are obtained 
if plumes are given a Gaussian profile in time, with randomly 
distributed peak times.) The c hoice of 15 minutes is consis- 
tent with the observations of IWelsch et al.l (1201 ih for flows 
on this spatial scale. 

Rotational flow. Here the velocity field v is a superposition 
of vortices 



i=l 



(X- 



■Xi)^ 
: v^i. 



(10) 



where (X/, Yi) is the centre of the vortex, r/ = 
^J(x- Xi)^ (y - Yi)^, and the vortex has an azimuthal ve- 
locity 

a. ( -rf . ^^^^ 



I ' 

Now 0) ^ and A = 0. We choose N = SI vortices, with 
centres (X/, ?/) on an extended grid as in model 1. Since 
oj and A are observed to have comparable scales (Sect. O, 
the widths Ri are again selected from a normal distribution 
with (Ri) = 1 Mm and sd(^/) = 0.3 Mm. The strengths 
are selected from a normal distribution with zero mean and 
sd(Q/) = 0.00048 s-^ The latter value makes the RMS oj 
in the box approximately equal to the RMS A in model 1 
(approximately 2 x 10"^ s"^). Again, a constant pattern of 
vortices is applied for each coherence time. 
3. Diverging and rotational flow. The velocity field is the su- 
perposition of those from models 1 and 2. 

In each model, an initial grid of 864 x 864 tracer points has been 
integrated in the same way as for the observational data, follow- 
ing trajectories of the model flows. The resulting a fields for 
each model after 6 hours (corresponding to 24 coherence times) 
are shown in Fig.[7^-c. The time evolution of J log^g Qdxdy is 
shown by the thin lines in Fig.[6]for models 1 (dotted), 2 (dashed) 
and 3 (solid). Note that the combined model 3 leads to a linear 
increase at a rate comparable to the observed flow (thick solid 
line). This rate is approximately double that of models 1 or 2. 



4.2. Origin of the Observed Pattern 

From the cr fields for the three analytical models in Fig. [71 we can 
see that the qualitative structure of the observed a field (Fig.O 
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(a) ti = 6 hours 
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Fig. 7. Finite-time Lyapunov exponent (FTLE) field cr after 6 hours for photospheric flows in the analytical models with (a) diverging 
flow only, (b) rotational flow only, and (c) both diverging and rotational flow. 



reflects predominantly the influence of divergence A in the ve- 
locity field. This is responsible for generating the visible ridge 
network of LCSs (or QSLs), interspersed with regions where 
(T < 0. Note that this network diff'ers from that of the magnetic 
flux distribution (Fig.[T|): the latter collects along the stable man- 
ifolds rather than the unstable manifolds. 

Model 2, with purely rotational flow (Fig. [TJ)), generates a 
very diflferent cr field: there are still localised ridges, but these are 
more diflfuse and space-filling than in the diverging flow. There 
are no regions of cr < 0. When both flows are superimposed 
(Fig. [7]:), the cr field visually resembles that of the diverging 
flow, with the contribution of vorticity being to locally perturb 
the ridge network and to infill some of the cr < regions with 
new ridges (although, in fact, J log^g Qdxdy increases at twice 
the rate of model 1). This model, where A and co are compara- 
ble in magnitude, corresponds most closely with the observed 
velocity field. 

Although rotational component of the combined flow has lit- 
tle visual eflfect on the LCS/QSL pattern, this does not mean that 
it is incapable of leading to gradients and subsequent reconnec- 
tion in the coronal magnetic field. Rather, the superposition of a 
diverging flow has perturbed the picture of the cr and Q fields so 
as to mask the contribution from oj. It may therefore prove use- 
ful for future analysis to extract the two components from the 
observed velocity field. This may be done through a Helmholtz 
decomposition 

V = Vdiv + Vrot, (12) 

where Vdiv = V0 and Vrot = V x (j/^e^), for scalar functions 
0(x, y, t), i//(x, y, t). This decomposition is not unique but we can 
fix a particular solution by specifying the boundary conditions 
that n • Vrot = on the boundary of our square region, where n 
is the unit normal to the boundary. Then v • n = Vdiv • n on the 
boundary. The function is determined by solving the Poisson 
equation = A with the Neumann boundary conditions that 
n • V0 = n • V on the boundary. The function il/ is determined by 
solving the Poisson equation 

VV = -co, (13) 

in this case with the Dirichlet condition of constant il/ on the 
boundary. 

Given co from the observed velocity, we calculated j/^(x, t) 
at each time sli ce by solving Eq uation ([T3]) using a standard fast 
Poisson solver dvan Loanlll992l) . The cr fields were then calcu- 
lated for Vdiv and Vrot separately, and are shown at 02:24UT in 



Fig. O This confirms our qualitative picture from the analyti- 
cal model of the rather diff'erent cr fields for the two types of 
flow. The thick dotted and dashed curves in Fig. [6] show that 
J logio Qdxdy for both Vdiv and Vrot increases at about half the 
rate of that for the combined velocity, again in accordance with 
the model. 

4.3. Predicted Structure from Faster Flows 

Since our method relies on the photospheric velocity field 
v(x, y, t) as input, the inferred magnetic structure is necessarily 
aflfected by limitations in the velocity observations used. In addi- 
tion to the limited spatial resolution of SOT (which has a diff'rac- 
tion limit of approximately 0.3"), the velocities we report have 
been averaged in both space and time in order to reduce noise in 
the local correlation tracking. Previous studies of granular con- 
vection on the Sun report cells of size 0.5-2 Mm with typical 
l ifetimes of ~ 10 minutes , and velocities from 0.5 - 1.5 km s"^ 
(Rieutord & Rinconl20 10'). While the "cells" we find are compa- 
rable to the larger end of this size range, our velocities (Fig.|2j)) 
are a factor of 5 - 10 smaller due to the eflfect of the averaging. 

To predict the magnetic structure that would arise from faster 
flows, we re-run the analytical model with plume velocity in- 
creased first to {Vi) - 0.3 km s"^ and then to (V/) = 0.5 km s"^. 
To keep the RMS vorticity comparable to the RMS divergence, 
we also increase sd(Q/) to 0.00144 s"^ and 0.0024 s"^ respec- 
tively. The resulting cr fields in models 1 , 2 and 3 after 6 hours 
are shown in Fig. [9l The results are striking: in model 1 (di- 
verging flow), the efl'ect of faster flows is to sharpen the LCSs, 
increasing the maximum value of cr but maintaining the overall 
pattern. By contrast, in model 2 (rotational flow), faster flows 
lead not only to sharper LCSs but also to greater filling of space 
with these structures. Interestingly, model 3 (combined flows) 
sees not only the sharpening efl'ect of model 1 but also an in- 
creased infilling of the cr < regions with new LCSs resulting 
from the vortices. 

The infilling arises because, with strong enough vorticity, 
LCSs "wind up" around centres in the flow pattern. This is par- 
ticularly evident in regions where the vorticity is strong but the 
divergence is weak; for example, compare the region around 
X - \,y = 1.5 in Figs.|9h-f- Th is "wind up" phenomenon was 
found by 'Pemouli n et alT (ll996) who computed QSLs in analyt- 
ical flux tubes of increasing twist. Similarly, Birn et al. (1989) 
found that field line connectivity in a 3D plasmoid varies on 
smaller and smaller scales as the axial field is reduced relative 
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Fig. 8. Finite-time Lyapunov exponent (FTLE) 
field cr at 02:24UT for the (a) diverging compo- 
nent Vdiv and (b) the rotational component Vrot 
of the Helmholtz-decomposed observed veloc- 
ity field. 



to the azimuthal field, with the boundaries between regions of 
diflTerent connectivity wrapping round and round the flux tube 
axis. 

The initial rate of increase in J log^g Qdxdy (correspond- 
ing to the slopes in Fig. [S]) increases from ~ 8 when (V/) = 
0.1 kms"^ to ~ 44 when (V,) = 0.3 kms"^ and ~ 91 when 
(Vi) = 0.5kms"\ although the rate of increase slows consid- 
erably in the latter case during the 6 hours. This may correspond 
to a "saturation" of the infilling evident in models 2 and 3 once 
the cr field has become homogenised, in accordance with previ- 
ous studies of LCSs in turbulent convection (ILapeyre 2002). At 
this stage, magnetic field lines in the whole region will have be- 
come mixed/braided in a manner likely to promote reconnection 
and subsequently heating of the coronal plasma. A saturation is 
also seen in the maximum value of cr, and this occurs after a 
shorter time for higher flow speed. This may be due to the LCS 
widths falling below the tracing grid scale, but it is unclear what 
other eflTects might cause such saturation; this bears further in- 
vestigation. 

Note that certain consequences of a faster plume velocity, 
namely (i) sharper LCS with higher peaks of cr, and (ii) a faster 
rate of increase in J log^Q Qdxdy, would also result if one left 
the plume velocities unchanged but increased their coherence 
time. This is demonstrated in Fig. [TOl which shows model re- 
sults after six hours with a single velocity pattern, rather than 
changing the pattern every 15 minutes. The peak cr (c.f. Fig. [7]) 
and slope of the J log^g Q dxdy curve (not shown) become com- 
parable to the run with (V/) = 0.3 km s"^ Yet there are fewer 
LCS, filling less of the area. This illustrates how the pattern of cr 
depends on the time history of the flow, not just on its pattern at 
any given instant. 

For simplicity, the models here are limited to convective 
cells with a single spatial scale, flow speed and lifetime. On 
the real Sun, convection seems to operate simultaneously on a 
range of scales. Experiments with superimposing a "supergran- 
ular" flow in the model (cells 10 times larger, with slower speeds 
and longer lifetimes) indicate that the LCS pattern and increase 
of J logiQ Q dxdy are determined primarily by the original, faster 
flow component. The supergranular flow can generate localised 
magnetic gradients only over a longer timescale of many hours. 



5. Conclusion 

We have proposed a practical method for inferring the topology 
of the 3D coronal magnetic field not by extrapolation (as is typ- 
ically used) but rather by integrating trajectories of an observed 
sequence of horizontal flows in the photosphere. Assuming an 



ideal evolution in the corona, this enables us to infer the re- 
sulting magnetic field line mapping between photospheric foot- 
points, and therefore the squashing factor Q of this mapping. 
The method has been demonstrated using a particular sequence 
of observed flows derived from Hinode SOT magnetograms, al- 
though it applies equally to any velocity data and does not re- 
quire magnetic information. Since the inferred field line map- 
ping results directly from the observed photospheric flow, quasi- 
separatrix layers in the coronal magnetic field correspond to so- 
called Lagrangian Coherent Structures in the flow. These are 
ridges of high finite-time Lyapunov exponent cr where neigh- 
bouring trajectories diverge most strongly. 

With a simple analytical model, we have demonstrated that 
the pattern of the cr (or equivalently Q) field differs significantly 
between a flow composed of irrotational convective plumes and 
one composed of purely incompressible vortices. The diverg- 
ing flow pattern in the first case leads to a network of long thin 
LCSs, whereas the vortical flow leads to a space-filling pattern of 
shorter LCSs that look a priori more favourable for widespread 
reconnection. The observed velocity field is found to have com- 
parable |V- v| and |Vxv|, and the observed cr field is in qualitative 
agreement with a combined model incorporating both eflTects. In 
the combined case, the appearance of the cr and Q fields follows 
that of the diverging flow model: the diverging part of the veloc- 
ity acts to quickly stretch and deform the picture. However, the 
rate of increase of integrated log^g Q in the combined model is 
double that of the original model, and we hypothesise that the 
vortical structure remains "hidden" in the magnetic field topol- 
ogy. We have demonstrated how the vortical part may be ex- 
tracted from an observed velocity field, but further study is re- 
quired to determine whether the vortical part is a more appropri- 
ate predictor of subsequent reconnection. 

Due to the limitations of the observational technique used for 
the demonstration in this paper, the typical flow speeds measured 
were only 0. 1 km s"\ a factor of 5 - 10 slower than real granular 
flows. From investigation of the analytical model, we predict that 
faster flow speeds (for the same size and lifetime of granules) 
will result in significantly faster development of strong gradients 
in the magnetic field. Our initial results are therefore very much 
a lower bound for the complexity that we expect to develop in the 
coronal magnetic field over this time. The model also indicates 
that, if the real vorticity is also larger, then the combined cr field 
will show greater infilling of LCSs. In this case, the mixing of 
trajectories is sufl&cient that the model begins to show a process 
of "homogenisation" of the cr field as found in simulations of tur- 
bulence ( Lapeyre 2002). A proper investigation of the rate of this 
mixing will require higher resolution velocity data, but the sim- 
ple model indicates that it is likely to be eflfective on timescales 
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Fig. 9. Effect of an increased flow speed in the analytical model. In panels (a)-(c) the mean plume velocity is 0.3 kms \ while in 
panels (d)-(f) it is 0.5 kms"^. As in Fig. [71 the cr field is shown after 6 hours for models 1 (left column), 2 (middle), and 3 (right). 
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Fig. 10. Finite-time Lyapunov exponent (FTLE) field cr after 6 hours of a stationary velocity field in the analytical models with (a) 
diverging flow only, (b) rotational flow only, and (c) both diverging and rotational flow. 



much shorter than the observed 12 hours. Realistic numerical 
simulations of photospheric convec tion (e.g.. Stein & Nordlund 
ll998l:lGudiksen & Nordlundll2005l) could also give tighter con- 
straints on the expected magnetic structure. 

Finally, the method proposed here — which assumes a per- 
fectly ideal evolution in the corona — will break down if and 
when sufl&ciently high magnetic gradients have formed for mag- 
netic reconnection to set in. Determining this threshold for re- 
connection will likely require detailed study of numerical MHD 
simulations. To this end, we have proposed a method for recon- 
structing a 3D magnetic field with the inferred field line map- 
ping. This field is neither unique nor (likely) in equilibrium, but 



it has the significant advantage over existing extrapolation tech- 
niques of having the correct field line topology. 
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to make B divergence-free and match a given normal component 
Bz(x, J, 0) on the lower boundary. From V • B = we find 



— = -M, 

dz 



(A.2) 



where the z-derivative is taken along a magnetic field line and 
A = V • V as before. Integrating from ^ = to ^ = z along the field 
line f(xo,jo,0 gives 



3) 



n/l ^1 nz 

— =- A(f(xo,yoJ))dt, (A. 

JB,(xo,yo,0) ^ Jo 

SO at the point (x, y, z) = f (xq, yo, z) we can integrate to find 

A(x, y, z) = B,(xo, Jo, 0) exp |- A(f (xq, Jo, O) dt^ . (A.4) 

Notice that the magnetic field is entirely determined by knowing 
both its magnetic field lines and the distribution of on the 
lower boundary. The magnetic field in Fig. [T] was constructed in 
this way from the observed velocity field. 



Appendix A: Magnetic Field Construction 

The method outlined in this paper infers the magnetic field line 
mapping from a photospheric flow. However, the field line map- 
ping defines the 3D magnetic field B only up to an ideal defor- 
mation. In particular, a hypothetical plasma flow in the volume, 
which vanishes on the photospheric boundaries, can deform B 
while leaving the field line mapping invariant. Conversely, two 
fields B resulting from applying the same photospheric footpoint 
motions to the same initial magnetic field can diff'er at most by 
an ideal deformation. 

It would be useful for future investigations to generate a par- 
ticular magnetic field of the required topology. Here we present 
a general method for constructing such a field. Given a velocity 
field v(x, y, t), the strategy is to set 

B(x,y,z) = A(x,y,z){vjc(x,y,z)e^ + Vy(x,y,z)ey e,) (A.l) 

so that the field lines of B are simply the trajectories of v with 
z corresponding to time. The scalar function A is then adjusted 
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